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"TRANSONIC  SWEPT-WING  ANALYSES  USING  ASYMPTOTIC 
AND  OTHER  NUMERICAL  METHODS" 

(AIAA  Paper  80-0342) 
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Page  1,  left  column,  17th  line  in  Abstract,  change  "swept-forward-wing"  to 
"swept- forward  wing". 

Page  1,  right  column,  15th  line  from  bottom,  insert  "and"  between  "wings" 
and  "conven-". 

Page  I,  right  column  4th  line  from  the  bottom,  change  "analysis"  to 
"analyses". 

Page  2,  left  column,  Eq.  (2.1),  change  "M^"  to  "M  ". 

Page  2,  left  column,  !4th  line  from  the  bottom,  replace  "that"  with  "the 
scale". 

" (57) "  " (56) " 

Page  2,  right  column,  end  of  second  paragraph,  change  '  to  '  ’  ; 

below  Eq.  (2.6b),  change  "where"  to  "while";  Eq.  (2.6c),  replace  "m" 
with  "m". 


Page  3,  left  and  right  columns.  Section  2.2,  Weissinger  has  been  cited  and 
misspelled  four  times. 

Page  3,  left  column,  5th  line  below(Eg.  (2.7K  add  comma  after  "turn"; 

7th  line  below  Eq .  (2.7),  replace  '*7/  by  '(36)".  second  line  of  &  2.2, 
move  "28-30"  to  the  front  of  the  parenthesis. 

C  iif  ■jQV1  "  ( " 

Page  3,  right  column,  13th  line  of  ?  2.3,  replace  with  'J/  . 

Page  4,  left  column,  24th  line  from  the  bottom,  insert  "will"  between  "but1 
and  "reduce". 

Page  6,  right  column,  B  3-4,  2nd  line  from  the  top,  change  "admits"  to 
"admit". 


Page  6,  right  column,  Eq.  (3-11),  replace  " "  w ‘  t  h  "  "• 

Page  7,  left  column,  t  3-5,  2nd  line  from  the  top,  change  "existance"  to 
"exi stence" . 


<L 


Page  8,  left  column,  23rd  line  from  the  top,  change  "because"  to  "from". 

M ( 53) H  M ( 52) " 

Page  8,  left  column,  footnote,  replace  with  p  '  . 


Errata  (continued) 


Page  9,  left  column,  §  4 ,  3rd  line  from  the  top,  change  "lead"  to  "leads". 

Page  9,  right  column,  4th  line  from  the  top,  change  "posess"  to  "possess". 

Page  11,  left  column,  7th  line  from  the  top,  change  "vort i c i t i es"  to 
"vort ici ty". 

Page  11,  left  column,  26th  line  from  the  too,  change  "30o"  to  "22.5°"  and 
"M  =  .755"  to  "M  =  .824". 

CD  CD 

Page  11,  left  column,  39th  line  from  the  top,  change  "Fut hernore"  to  "Furthermore 

Page  il,  right  column,  3rd  line  from  the  bottom,  change  "convergenced"  to 

"converged". 

J'lfcnV  Mfcqpi 

5.2,  replace  '  with  v ’  . 

Pane  12,  left  column,  2nd  line  of  §  5-3  replace  "(4.62)"  with  "(11,62)"; 

3rd  paragraph  .§5.2,  4th  line,  add  Ref.  22. 

Page  13,  left  column,  2nd  line  from  the  end  of  §  5,  change  "is"  to  "based". 

Page  13,  loft  column,  30th  line  from  the  top,  delete  "is". 

Page  13,  left  column,  22nd  line  from  the  bottom,  change  "vor t i c i t ies"  to 

"vort i c i ty". 

Page  13,  right  column,  5th  line  from  the  bottom,  correct  the  author's  name 
to  read  "Kuchemann". 

Page  13,  right  column,  Reference  1:  change  "Busenan"  to  "Buserann";  write 
"Uberscha  I  1  ga-.' ca  Windigkeit"  as  a  single  word. 

Page  13,  right  column.  Reference  4:  change  "Kucheman"  to  "Kuchemann". 

Page  14,  left  column.  Reft  rence  14a:  change  "to  appear  in...."  to  "Journal 
of  Fluid  Mechanics','  Vol.  97,  pt.  3,  1980,  pp.  531_555", 

°age  14,  left  column.  Reference  19:  write  "Trafluqel  Theorie"  as  a  single 
word . 

Page  14,  right  column,  Reference  28:  change  "Wei  singer"  to  "We i ss i nger" . 

Page  15,  right  column.  Reference  55:  delete  "and"  from  the  3rd  line. 

Page  15,  right  column,  Reference  57:  iri3ke  the  first  letters  in  "Senkrechten", 

"Einer",  and  "Gekrummten"  lower  case. 

Page  15,  right  column,  the  Caption  of  Fig.  1,  3rd  line  from  the  bottom,  change 
"swept-wing"  to  "swept-wings". 

Page  16,  the  Caption  of  Fig.  7.  change  "D"  to  lower  "d". 


Errata:  (1980),  vol.  97,  pt.  3,  pp.  531-556,  June  4,  1980 

Cheng  and  Meng. 

Page  531,  third  line  of  the  "Introduction”:  add  "in  a  transonic 
or  supersonic  stream"  after  "airfoil  section". 

Page  535,  the  line  preceding  Eq.  (2.7) :  Replace  "y"  in  a  and  b  by  "y" 

Page  545,  second  line  of  footnote:  Replace  "tfx"  by  "¥?c" . 

Page  546,  Eg.  (4.14b),  second  line:  Change  |<*j2  to  }fj4. 

Page  549,  2nd  and  3rd  lines  be lew  Fig.  3:  Change  "induced"  to  "anti- 
symmetric  part  of  the  induced";  delete  "(CL<0)"  and  "(Cj>0) 

5th  line  from  the  bottom:  Replace  by  "$>2  • 

”1  A  "  "1  rJ  " 

Figs.  4(a-c)  and  Fig  5:  change  ^(x+l)  to  ^-(x+1) 
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Transonic  swept  wings  are  anal)  red  as  a  lifting-line  problem  under  a  small-disturbance  approximation.  Basic 
concepts  and  principal  results  of  the  asymptotic  theory  are  discussed.  The  study  focuses  on  straight  oblique 
wings  and  V-shaped  swept  wings,  of  which  the  local  centerline  curvature  can  be  equated  to  zero.  The  three- 
dimensional  1 3-D)  perturbation  of  the  nonlinear  component  flow  admits  a  similarity  flow  structure  but  requires 
that  all  wing  sections  are  generated  from  a  single  airfoil  profile;  the  reduced  2-D  problems  in  this  case  are  solved 
only  once  for  all  span  stations.  Examples  of  solutions  involving  high  subcritical  and  slightly  supercritical 
component  flows  are  demonstrated  and  compared  with  surface  pressure  data  from  3-1)  computer  codes  based  on 
the  full-potential  equation  (HO  22).  Except  in  the  neighborhood  of  leading  edges,  where  the  small-disturbance 
assumption  breaks  down,  and  in  the  vicinities  of  wing  lips  and  the  symmetry  plane,  where  neither  the  theory  nor 
the  3-D  codes  may  claim  full  validity  ,  reasonable  agreement  is  consistently  found.  The  explicit  results  from  the 
upwash  analysis,  along  with  the  similarity  flow  structure,  provides  a  rational  approach  to  the  control  of  3-D 
effects  in  transonic  aerodynamic  design  studies. 


I.  Introduction 

UR  understanding  in  the  aerodynamics  of  wing  sweep 
and  its  use  to  control  the  compressibility  effect  have  been 
limited  in  the  past  mostly  to  problems  in  the  linear  flow 
regimes. 15  Current  interest  in  aircraft  wing  design  has 
focused  on  adopting  two-dimensional  (2-D)  supercritical 
airfoil  data,6'11  where  the  flow  in  this  domain  is  nonlinear 
with  a  transonic  flow  component.  The  problems  of  sweep  and 
the  3-D  corrections  to  mixed  component  flows  have  recently 
been  studied  by  Cheng,  Cook,  and  others 12  20  as  extensions  of 
the  classical  lifting-line  theory. 21  22  As  a  sequel  to  Cheng  and 
Meng’s  work, 1,16  this  paper  will  present  a  more  complete 
discussion  of  its  applications,  considering  examples  nol 
previously  analyzed.  An  essential  part  of  th  study  is  com¬ 
parisons  with  solutions  generated  from  ex>  .mg  full-potential 
computer  codes;  certain  apsects  for  fu'  jie  improvement  of 
the  latter  will  also  be  brought  out . 

In  the  domain  of  interest,  the  tiowfield  far  from  the  wing 
section  should  pertain  to  a  high  subsonic,  or  linear  sonic, 
outer  flow,  which  is  representable  in  the  leading  ap¬ 
proximation  by  a  solution  to  the  Prandtl-Glauert  equation 
corresponding  to  a  swept  lifting  line.  Thus,  the  theoretical 
treatments 12  20  mentioned  have  genesis  in  Prandtl’s  original 
lifting-line  idea,2122  even  though  the  corresponding  inner 
problem  is  nonlinear  and  may  involve  an  embedded  super¬ 
sonic  region.  Apart  from  gaining  a  better  physical  insight  and 
a  greater  simplicity  in  the  analysis,  the  approach  reflects  a 
desire  for  implementing  the  current  computer -oriented  3-D 
methods21 26  which,  though  being  very  powerful,  share  the 
well-known  problem  of  trading  accuracy  with  computer  cost 

Received  Aug.  20.  19X0;  revision  received  March  IX.  1981. 
Copyright  ©  American  Institute  of  Aeronautics  and  Astronautics, 
Inc..  1981 ,  All  rights  reserved. 

•Professor,  Dept  of  Aerospace  Engineering.  Member  AIAA. 

♦  Research  Assistant;  presently  at  Hydrodynamics  Division. 
Rockwell  International  Corporation. 

{Staff  Scientist,  Research  Department.  Member  AIAA. 

^Senior  Research  Scientist,  Aeronautics  Division.  Member  AIAA. 


and  storage  requirements.  One  obvious  contribution  to  which 
the  approach  may  lead  is  providing  an  asymptotic  basis  upon 
which  numerical  results  from  current  3-D  codes  can  be 
correlated  and  their  adequacy  tested. 

The  asymptotic  theory  underlying  the  present  study  has 
been  given  quite  thoroughly  in  Cheng  and  Meng's  work. 11 16 
A  similar  formulation  has  been  developed  by  Cook  for 
straight  oblique  wings. 211  Most  high  aspect  ratio  wing  analyses 
in  classical  literature  arc  concerned  with  unyawed,  straight 
wings,  and  were  tied  closely  to  Prandtl’s  lifting-line  theory: 
documented  works  in  this  category  are  numerous  (see  review 
by  Jones  and  Cohen1).  Some  of  the  earlier  works10  ”  possess 
features  which  are  quite  relevant  to  the  present  study,  and 
have  been  reviewed  along  with  transonic  lifting-line  works  r  20 
in  Ref.  34. 

Section  II  will  discuss  the  basic  concepts  of  the  theory  and 
their  physical  implications.  The  distinction  from  classical 
lifting-line  analyses  will  be  apparent.  In  order  to  allow  a  fuller 
discussion  of  the  solution  examples  and  their  comparison  with 
full-potential  solutions  (see  Sec.  IV),  mathematical  details 
described  in  Sec.  Ill  will  be  kept  to  a  minimum,  of  which 
readers  are  referred  to  the  published  works. 14 16 

A  limitation  of  the  present  approach  is  the  assumption  of  a 
smooth  centerline,  which  obviously  breaks  down  in  the 
vicinity  of  the  apex  of  a  conventional  (symmetric)  swept  wing. 

The  example  analyzed  will  show,  however,  that  this  region  of 
nonuniformity  is  limited  in  extent;  agreement  with  the  full- 
potential  solution  remains  reasonable  at  span  stations  well 
within  one  root  chord  from  the  symmetry  plane.  An  im¬ 
portant  property  uncovered  from  the  study  of  Refs.  14-16  is 
the  similarity  flow  structure  which  permits  the  reduced 
problem  to  be  solved  only  once  for  all  span  stations.  In¬ 
terestingly,  this  properly  remains  even  after  the  small- 
disturbance  assumption  is  removed,  as  demonstrated  by  a 
formulation  based  on  full-potential  equations  in  Ref.  36. 

The  following  analysis  will  deal  mainly  with  centerlines 
composed  of  one  or  two  straight  segments  such  as  a  straight 
oblique  wing  or  a  V-shaped  swept  wing,  respectively.  All 
examples  studied  'B*1  *■  1  **■  ** — *  * 
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Fig.  I  C  artesian  and  orthogonal  curvilinear  coordinates  in  the  wing 
plane  ;  =  =0. 


II.  Essential  Elements  of  the  Theory 

Following  Prandtl's  lifting-line  theory  for  high  aspect  ratio 
wings,’1  two  (asymptotically)  distinct  flow  regions  are 
considered:  1)  a  nearly  planar  (2-D)  region  next  to  the  wing 
section  with  a  streamwise  length  scale  comparable  to  the 
typical  wing  chord  c„;  2)  a  fully  3-D  domain  with  its  lateral 
si/e  comparable  to  the  wing  half-span  b.  These  are  the  inner 
and  outer  regions  of  the  singular  perturbation  theory.”  Since 
the  aspect  ratio  M,  as2b/c„  is  high,  the  wing,  along  with  its 
near  wake,  is  perceived  in  the  outer  flow  region  as  a  lifting 
line — a  line  of  singularities,  to  be  more  precise. 

Students  of  classical  aerodynamics  are  familiar  with  the 
idea  that  the  velocity  induced  by  the  trailing  vortex  sheet 
causes  a  change  in  the  local  flow  angle,  proportional  to  /R;  ' . 
Implicit  arc  the  assumptions  that  the  flow  around  the  wing 
section  is  strictly  planar,  and  that  the  induced  velocities 
defined  in  the  classical  sense  remain  bounded.  These 
stipulations  do  not  hold  for  a  swept  wing  with  a  nonvanishing 
sweep  angle.  The  distinct  features  brought  out  in  the  analysis 
for  swept  wings  are  described  below . 

Centerline  Curvature  and  Spanwise  Density  Variation 

The  partial-differential  equation  governing  the  component 
flow  around  the  wing  section  must  be  corrected  for  com¬ 
pressibility  effects  and  for  the  3-D  effect  due  to  a  non¬ 
vanishing  local  curvature  of  the  centerline  of  the  planform, 
dA/dy'  (see  Fig.  1  for  coordinates  and  definition  of  A).  The 
former  includes  a  correction  resulting  from  the  spanwise 
density  variation  (see  Sec.  III). 

Spanwise  Component  of  Wake  Vorticity 

Owing  to  the  sweep,  the  wake  vorticity  just  behind  the 
trailing  edge  has  a  nonvanishing  component  along  the  cen¬ 
terline.  Thus,  the  tangential  velocity  component  normal  to  the 
centerline  has  a  nonvanishing  jump  across  the  trailing  vortex 
(TV)  sheet.121’  Therefore,  the  solution  to  the  reduced  inner 
problem  is  no  longer  planar,  but  must  be  corrected  in  both  the 
partial-differential  equation  and  in  the  boundary  condition  on 
the  wing  trace,  i.e.,  the  TV  sheet. 

I.ngarilhmic  Cpwash 

The  induced  upwash  effect  on  the  wing  section  charac¬ 
teristics  is  determined  from  the  behavior  of  the  outer  (lifting¬ 
line)  solution  in  the  vicinity  of  the  lifting  line.  In  approaching 
the  latter,  the  outer  solution  is  dominated  by  that  of  a  line 
vortex  representing  the  effect  of  bound  vortices  on  the  wing. 
Unlike  the  classical  theory  for  unyawed  straight  wings,-’1” 
the  induced  upwash  (obtained  after  subtracting  out  the  vortex 
singularity)  approaches  infinity  like  the  logarithm  of  the 
distance  from  the  lifting  line.  This  singular  behavior  is  needed 
for  matching  with  the  corresponding  (far-field)  behavior  of 


the  inner  solution  [see  Eq.  (14b)).  The  strength  of  this 
singularity  is  proportional  to  tanA(dI7dy ' ) ,  where  T  is  the 
circulation,  hence,  to  the  sweep  angle  and  the  shedded  vor- 
licity  at  the  span  station  in  question.  Upon  matching  with  the 
inner  solution,  the  logarithm  yields  an  additional  term  to  the 
upwash  proportional  to  ln(2b/c„)  =(n  which  dominates 
over  the  remaining  finite  part  of  the  upwash,  tV, .  In  view  of 
its  linear  dependence  on  df/dy ' ,  it  may  rightfully  be  regarded 
as  a  contribution  of  the  near  wake. 

The  origin  of  the  logarithmic  singularity  can  be  understood 
from  a  decomposition  of  the  locally  shedded  trailing  vorticity 
into  two  components,  parallel  and  normal  to  the  centerline 
(.v'=0).’4  The  component  parallel  to  x'  =  0  gives  a 
logarithmically  infinite  upwash  as.v'  —  0  *  . 

Logarithmically  large  upwash  corrections  occur  also  in 
unsteady  problems  and  in  problems  involving  curved  cen¬ 
terlines.  Analyses  of  the  outer  flow  behavior  near  the  cen¬ 
terline  has  been  described  for  a  straight  oblique  wing  in  Ref. 
12  for  incompressible  flow,  and  for  a  wing  with  straight  as 
well  as  curved  centerline  in  the  transonic  regime  in  Refs.  15 
and  16.  (Results  obtained  from  the  curved  lifting-line  theory 
for  incompressible  flows  in  the  steady  and  unsteady  cases 
have  not  yet  been  published.) 

Transonic  .Swept- Wing  Problems 

For  problems  involving  a  transonic  flow  component 
(Af„  =  A/„cosA~  I),  the  nonlinearity  familiar  in  the  transonic 
flow  study  must  be  adequately  treated  in  the  equation 
governing  the  inner  flow.  If  the  freestream  is  also  transonic, 
i.e.,  M„~  I,  the  conditions  A/„  - -  1  implies  that  the 
swept  angle  must  remain  sufficiently  small  in  the  nonlinear 
regime.  Denote  the  absolute  angle  of  attack  or  the  thickness 
ratio  of  the  wing  section,  whichever  is  larger,  by  a.  The  range 
of  sweep  must  be  such  that  (I  -  Af ; ) /or  '  =0(1),  in  order  to 
keep  the  component  flow  nonlinear.  This  requires 
A  =  0(a;  ').  Interestingly,  the  smallness  in  A  renders  the 
effect  of  the  spanwise  component  of  the  trailing  vorticity 
negligible  in  the  off-wing  boundary  condition,  but  the 
logarithmic  upwash  associated  with  this  spanwise  trailing 
vorticity  remains  a  dominating  feature  in  the  transonic 
problem. 

III.  Reduced  Transonic  Problems  and  Key  Equations 

Basic  Parameters  and  (ioveming  Equations 

We  consider  high  subsonic  flight  near  the  speed  of  sound. 
The  component  flow  in  this  case  is  necessarily  nonlinear  and 
may  become  supercritical.  If  the  wing  sections  are  assumed  to 
be  thin,  this  component  flow  at  each  span  station  will  be 
controlled  by  the  transonic  similarity  parameter  based  on  the 

component  Mach  number  A/„  mA/„cosA 

*„•(/  -Mi)/a}  'M„  (1) 

familiar  from  the  transonic  small-disturbance  theory.  ’’  The 

sweep  range  of  interest  is  that  which  will  keep  A'„  at  order 
unity,  while  A/„  does  not  exceed  I .  This  requires 

0-  ~7=0(/)  (2) 

a ' 

The  parameter  0  may  be  referred  to  as  a  reduced  sweep  angle. 
Essential  to  the  development  is  a  reduced  aspect  ratio,  the 
reciprocal  of  which  is 


The  asymptotic  analyses  of  Refs.  14-16  pertain  to  the  limit 
r— 0.  holding  K„  and  0  (as  well  as  a'  'f  ')  fixed.  With  the 
additional  requirement,  <  I,  the  foregoing  requirements 
on  K„  and  G  can  be  replaced  by 
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e’<*„=0(/)  (4) 

We  shall  summarize  the  key  equations  for  the  reduced  inner 
problem  using  a  system  of  curvilinear  orthogonal  coordinates 
(x'.v'.z')  shown  in  Fig.  1.  Introduce  the  inner  reduced 
variables 

xm2x'/c„  /C0 

<-> 

ymy'/b  0  *20/ay'  U„c„ 


Perturbation  Analysis 

The  3-D  corrections  in  question  represent  a  small  petur- 
bation  of  a  2-D  transonic  component  flow;  formally,  we  write 

4>=(t>0(x ,z;y)  +fd>,  (x,z;y)  +  ...  (10) 

with  0U  and  0,  satisfying 

(l^-(Y+/)0()1]  ~  +  ~]0„  =  0  (11a) 


as  suggested  by  the  similitude  of  the  transonic  small- 
disturbance  theory,  with  c„/ 2  taken  to  be  the  reference  length 
scale  and  U„»UJ„cosA  as  the  reference  velocity. ,1M  In  the 
above  equation,  0  is  the  perturbation  velocity  potential  and 
a,  is  simply  a/M,. 

For  swept  wings  with  straight  centerlines,  the  full-potential 
equation  can  be  reduced  for  small  t  under  Eq.  (4)  to 14  16 

9  r  .  y  +  /  .  ,1  . 

2  ^  'J  +  =  (6) 


r  d1  d2  .  Ji. 

I  (*„ -iy+i)4>u, 1  rn  +  tt-,  -  (7+  0<<>, /u  —  U,  =2&<t>oii 

(  dz-  dx)  (Ilb) 

To  be  sure,  0„  is  governed  by  the  2-D  von  Kdrman  equation  in 
.V  and  z,  and  0;  satisfies  a  linear  nonhomogeneous  partial- 
differential  equation  (PDE)  in  x  and  i.  Throughout  Eqs.  (11), 
the  spanwise  variable  y  enters  merely  as  a  parameter.  The 
wing  boundary  conditions  are 


with  a  remainder  of  order  e- .  The  left-hand  members  are  the 
familiar  terms  in  the  2-D  von  Kdrmdn  equation,  r  and  the 
term  of  order  t  on  the  right-hand  side  is  the  3-D  com¬ 
pressibility  correction  for  the  spanwise  density  variation.  The 
pressure  coefficient  can  be  computed  from  0  as 


/90„  \  9  - 

(«)„“«**  <12a) 

/90, \  d  - 

(12b) 


--7“,  = -A'OS- A cry' 4,  O) 

with  a  remainder  comparable  to  Note  that  as  long  as 
O  =  0( I ),  C.  is  unaffected  directly  by  the  spanwise  velocity  0, . 
The  inclusion  of  M„  in  the  definition  of  A'„  leads  to  an  im¬ 
provement  in  predicting  the  critical  speed,  hence,  the  wave 
pattern  near  the  sonic  boundary,  as  ascertained  analytically 
by  Cheng  and  Meng. 1,1  We  point  out  that  most  cosA  factors 
appearing  in  Eqs.  (1-7),  except  that  in  (1-M2)=I- 
A/icos-’A  of  K„,  could  be  replaced  by  unity,  consistent  with 
the  asymptotic  theory.  However,  their  retention  will  allow  a 
greater  range  of  sweep  without  sacrificing  accuracy l?''6  (see 
also  Refs.  23  and  38-40). 

Inner  Boundary  Conditions 

A  nearly  planar  wing  will  be  considered,  of  which  the 
ordinate  is  expressible  as 

z'  ~C 2  («•?*  (x.y)  +a'  'Zgiy)  +  a<(x-x„)/(y)  |  (8) 

where  Z  *  ,  Z„,  and  /  are  assumed  to  be  of  unit  order  and  the 
"±”  signifies  the  top  and  bottom  surfaces.  The  unit-order 
functions  ZH  and  /  represent  an  upward  wing  bend  and  a  wing 
twist  sufficient  for  controlling  the  unwarranted  3-D  effects, 
such  as  the  asymmetrical  upwash  induced  by  the  wake  vor- 
ticity  behind  an  oblique  wing.'  The  impermeable  wing 
boundary  condition  can  be  transferred  to  the  wing  plane 
(c  =0),  subject  to  an  error  of  the  order  r\  as 


> 


(9) 


The  2-D  system  is  elliptic  or  hyperbolic,  depending  on 
whether  A„  -  (7  +  i)0«,  is  positive  or  negative.  The  coef¬ 
ficients  0,  and  .({’  must  allow  a  weak  dependence  on  <  to 
accommodate  the  logarithmic  upwash.  Where  a  shock  occurs, 
the  shock  boundary  is  described  as 

x =.vj’(c,y)  +  t.v '/(:;y)  + ...  (13) 

The  jump  conditions  transferred  to  the  unperturbed  shock 
boundary,  x=xl„>,  are  detailed  in  Refs.  14-16. 

Solution  Behavior  Far  from  the  Wing  Section  and  Matching 

Far  from  the  wing  section,  the  inner  equations  (II)  admit 
the  asymptotic  behavior  for  large  ^mx  +  isljt^z 

'(vfc-)+}sH+-  ,,4a) 

9  ,df„  .  . 

"4M 

The  remainder  in  Eq.  (14a)  is  dominated  by  terms  comparable 
to  a  doublet,  including  r„xf«  Ifl/  1  fl 2  and  f  !x! /  I  fl 4  (see 
Ref.  16,  Appendix  11).  The  remainder  in  Eq.  (I4bJ  is  com¬ 
parable  to  unity,  including  those  of  orders  (f«lfl)2  and 
(61  if  I ),  one  of  which  signifies  a  line  source. 1516 

The  potential  calculated  from  0o  +  (0,  based  on  Eqs.  (14), 
identified  as  the  outer  expansion  of  the  inner  solution,  can  be 
matched  with  the_  inner  expansion  of  the  outer  solution 
calculated  from  0„  mentioned  in  Sec.  II.  The  matching 
identifies  the  local  circulation  f„(y)  with  that  in  the  outer 
(lilting-line)  solution,  and  determines  the  upwash  correction 
C',  (  v)  from  a  finite  part  of  the  induced  velocity  at  the  lifting 
line  M',  as  "■'* 


where  the  subscript  w  refers  to  the  wing  surface. 

The  requirement  of  continuity  of  the  normal  velocity 
component  and  the  pressure  across  the  trailing  vortex  sheet 
(of  zero  thickness),  leads  to  the  condition  that  0,  and  0.,  not 
0 ,, ,  must  be  continuous  across  the  wing  trace.  Flcncc,  unlike 
the  swept-wing  problem  in  other  flight  Mach  number 
ranges, 12 11  the  inner  transonic  component  flow  can  be  treated 
without  the  vortical  wake,  subject  to  errors  of  the  order  t!. 


(15) 


The  inner  problem  is  thus  reduced  to  solving  Eqs.  (I  la)  and 
(lib)  with  boundary  conditions,  Eq.  (12),  and  the  far-field 
behavior,  Eqs.  (14),  using  the  results  of  matching,  Eq.  (15). 
The  Kutta  condition  is  essential  and  must  be  observed  at  the 
trailing  edge  for  0„  as  well  as  0;. 
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Fig.  2  Ip  wash  function  £  illustrated  for  \  1  -  V;  -Hf  =  8.39  at  three 
yaw  angles  for  an  elliptic  and  an  extended  span  distribution  in  t0. 


Induced  I'pwash  on  Swept  Wings — F.xamplesand  Nonuniformity 

To  provide  a  more  concrete  idea  on  the  induced  upwash 
distribution  and  its  difference  from  the  classical  result,  we 
shall  illustrate  the  function  IF',  [see  Eq.  (15))  for  straight 
oblique  wings  with  two  types  of  loads,  and  for  the  V-shaped 
swept-forward  and  swept-back  wings  with  elliptic  loads.  The 
function  IF',(v)  can  generally  be  written  in  two  parts 

W,  =  '/:  (£  +  £',)  (16) 

where,  with  mm  tan  A  =  tanA/VI  -  A/-,. 

4t£(.v)  =  -sinAr,',(  v)  [2MVA„/ 1)  +  2 


+  <’«((/- v-’)sec-’A(  |  +  f,;(.v)& 


I  —y  I  +  sinA 
/  +  v  /  -  sin A 


+  I  [  /  -  sinA(.v)  sgnfv,  -v)  )d.v,  (17a) 

■’  '  yt  -y 

and  t)  (v)  is  identically  zero  for  a  straight  oblique  wing, 
w  hereas  for  a  V-shaped  swept  wing,  it  can  be  calculated  from 


4  t£‘,  (y ) 


1  )  r  _ mj  l.v,  I  —  l.v  I  ) _ 

i  (v,  -  y )  •  L  v  m  -  (  I  v ,  I  -  Ivl )  •  +  ( v;  -  v)  • 


Fig.  3  I'pwash  function  £  +  £  ^  illustrated  for  elliptically  loaded 
swept-forward  and  swept-back  wings  of  At,  «  14  at  =0.755. 
Refer  to  text  for  conversion  In  other  aspect  ratio  and  Mach  numbers. 


/R,  is  9.18,  corresponding  to  Mm  =0.755  and  *  14;  results 
for  five  values  of  A  (0, ±22.5,  and  ±45  deg)  are  shown.  At 
=0.775,  A=±22.5  and  ±45  deg  corresponds  to 
A  =  ±  15.2  and  ±33.33  deg,  respectively.  The  need  for  a 
washout  on  a  swept-back  wing  is  quite  evident  from  the  large 
reduction  in  l£  +  £',l  near  the  tip  shown.  The  usefulness  of 
adopting  a  forward  sweep  ( A < 0)  to  revert  the  above  effect 
(and  to  enhance  maneuverability  in  roll)  is  also  quite  apparent 
from  Fig.  3.  In  the  outer  portion  of  the  swept-back  wing,  the 
upwash  behavior  is  similar  to  that  on  the  downstream  panel  of 
a  straight  oblique  wing.  Near  the  apex  (.1=0)  of  the  swept- 
back  wing,  the  induced  upwash  is  negative,  i.e.,  a  downwash, 
and  becomes  unbound  as  v— 0  like 

t)  ~ -£„(()) (sinA/dir.v)  (18) 


-sinA(.v)sgn(v,  -v)  Id.v,  (17b) 

where  the  subscript  0  refers  to  the  leading  (zeroth)  ap¬ 
proximation,  and  the  prime  signifies  a  derivative. 

The  asymmetrical  upwash  proportional  to  f,;(.v)  in  £  is 
responsible  for  the  unbalanced  rolling  moment  of  an  oblique 
(pivoted)  wing,  unless  twist,  wing  bend,  or  special  pivot 
location  is  introduced  in  the  design. ' 11  "  The  extent  to  which 
this  asymmetry  depends  on  the  pivoting  angle  A  and  on  the 
span  loading  £„(>•)  is  illustrated  in  Fig.  2,  where  the  upwash 
functions  £(y)  for  oblique  wings  with  an  elliptic  and  an 
"extended-span  distribution"  are  shown  for  A  =  0,  22.5,  and 
45  deg.  (The  extended  span  load  considered  has  a  root¬ 
bending  moment  equal  to  that  of  an  elliptic  load  for  the  same 
lift  but  lesser  drag41;  the  span,  however,  is  1.15  times  that  of 
the  elliptic  one.)  These  £\  were  computed  front  Eq.  ( 1 7a ) 
with  £(y)  normalized  by  the  midspan  value  of  £„,  and  the 
(reduced)  aspect  ratio  taken  to  be  such  that  VA„/f*s/l  -  M:„ 
/R(  »s/l  -  ;W^Af(sccA  »  8.4. 

The  induced  flow  angle  will  also  depend  on  the  centerline 
geometry.  Figure  3  shows  the  spanwise  distribution  of  the 
normalized  £  +  £)  for  symmetric  swept-forward  and  swept- 
back  wings  at  different  degrees  of  sweep.  It  is  assumed  in  the 
calculations  that  the  span  loading  in  the  leading  ap¬ 
proximation,  f„<» ,  is  elliptic,  and  that  the  product  v^l  -  M:m 


which  corresponds  to  the  velocity  induced  by  the  bound 
vortex  on  the  opposite  wing  panel.  (Compare  the  thin  dash 
curve  with  the  corresponding  £  +  £•,  in  Fig.  3.) 

Similarity  Flow  Structure 

Owing  to  the  linearity,  the  3-D  correction  can  be 
decomposed  into  separate  parts.  There  is  an  important  class 
of  wing  surface  geometry,  for  which  each  of  these  (properly 
scaled)  separate  parts  have  similarity  solutions  independent  of 
v.  as  does  the  basic  solution  4>„.  The  reduced  2-D  equation 
system  in  this  case  can  be  solved  once  for  all  span  stations. 
This  wing  class  requires  that  the  wing  section  at  each  span 
station  be  generated  from  the  same  airfoil  profile,  keeping  the 
same  thickness  ratio  and  incidence.  That  is,  the  function  of 
7.  *  (  v, y)  of  F.q.  (8)  is  required  to  conform  to 

7‘  =c(.v)/‘ <a7c)  (19) 

where  c(y)  ■c(.v)/c„.  Implicit  in  Eq.  (19)  is  that  the  cen¬ 
terline  v=;  =  0  is  the  common  straight  axis  for  the  similar 
wing  sections  at  different  span  stations.  The  location  of  this 
straight  axis  in  percentage  chord  may,  however,  be  arbitrarily 
set. 

The  similarity  variables  used  are 

Xmx/c  Zmz/C  CWCOsA.V  •>>’  4>m&/c 


(20) 
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Fig.  4  Spanwise  distribution  of  sectional  lift  coefficient  computed 
for  an  oblique  wing  by  solutions  to  the  full-potential  equation  (FLO 
22)  from  two  sources. 

The  similarity  flow  structure  admissible  under  Eq.  (19),  with 
due  allowance  for  the  upwash  correction  for  an  arbitrary 
loading,  is  given  by 

i  =  in(x.z)  +  eGc'0,(.?,c)  -f/K't'iZ 

+  t [ >/*''£', (v)  +  /(j it)  +  QZ'„ (y)  ] 0, (.v,f)  (21 ) 

where  c'  *dc/dy.  Z'H^dZB/dv.  There  is  a  similar  expression 
for  the  perturbed  shock  boundary.  I'"U4  As  indicated  earlier, 
0,i,  0,,  and  0,  are  independent  of  y.  Note  that  the  third  term 
on  the  right-hand  side  of  Eq.  (21)  is  introduced  so  that  the 
upwash  effect  is  transferred  through  the  solution  0;  to  the 
wing  surface  as  part  of  the  incidence  correction  to  0,  whereas 
the  second  term,  eG  r'0(,  accounts  only  for  the 
nonhomogeneous  correction  to  the  PDE  owing  to  the 
spanwise  density  variation. 

The  equations  and  boundary  conditions  in  the  reduced 
problems  for  the  three  similarity  solutions  0„,  0,,  and  0,  can 
be  obtained  readily  from  Eqs.  (II),  and  have  been  described 
in  detail  in  Refs.  14-16  and  34.  We  shall  not  repeat  it  here, 
except  to  point  out  that  the  z  derivatives  of  0,  and  0,  on  the 
wing  are  required  to  be  zero  and  unity,  respectively,  and  that 
the  upwash  function  C',  in  Eq.  (21 )  is  determined  from  C)  as 

0  f„  .  . 

C',  =  C',~  —  -fr'fcc  (22) 

2K  n  jr 

IV.  Examples  and  Comparisons 
with  Full-Potential  Solutions 

Inasmuch  as  the  existence  and  uniqueness  of  the  solutions 
cannot  be  easily  investigated,  demonstration  of  numerical 
solutions  to  the  reduced  problem  is  an  essential  part  of  the 
study. 

Remarks  on  Existing  3-H  Potential  Computer  Code 

The  crude  description  of  the  far  field  in  most  discretized 
flowfield  computation  methods  is  well  known.  This  problem 
becomes  more  serious  in  the  3-D  cases. 22  29  It  is  not  all  clear 
from  the  published  data  whether  the  grid  distributions 
used  therein  are  sufficiently  refined  for  the  purpose  of 
adequately  describing  the  upwash  induced  by  the  lar-wake 
vort icily,  which  is  crucial  in  the  analysis  of  high  aspect  ratio 
wings. 

There  appears  to  be  an  additional  problem  brought  about 
by  the  scarcity  of  the  span  stations  available  in  current  3-D 
potential-flow  programs. 2<  29  In  the  computer  code  lor  planar 


wings  based  on  Ref.  26,  referred  to  by  most  users  as  “FLO 
22,”  a  total  of  21  span  stations  are  allowed.  Applying  it  to  an 
oblique  wing,  for  example,  there  remains  only  10  stations  on 
each  wing  panel,  with  the  wing  tip  represented  by  a  trapezoid. 
Therefore,  one  cannot  attach  too  much  confidence  to  the 
numerical  data  obtained  for  the  tip  region. 1  It  has  also  been 
known  among  users  that  the  span  load  and  sectional  lift 
coefficient  are  far  from  being  smooth — the  values  at  suc¬ 
cessive  stations  appear  to  alternate  about  some  mean-value 
curve.  Figure  4  shows  typical  results  of  a  sectional  lift 
coefficient  obtained  from  the  converged  solutions  generated 
by  two  versions  of  the  FLO  22  code  for  a  6%  thick  elliptic 
wing  pivoted  at  22.5  deg  at  =0.824.  The  irregularity 
mentioned  is  clearly  noticeable  and  requires  some  caution  in 
making  comparisons.  The  code  is,  nevertheless,  known  to  be 
capable  of  reproducing  the  span  loads  determined  by  the 
(linear)  panel  methods  in  the  subsonic  range. 29 

Presumably,  this  irregularity  has  not  been  observed  in  the 
3-D  computations  using  the  transonic  small-disturbance 
(TSD)  equations24;  however,  existing  3-D  TSD  codes  are  not 
directly  applicable  to  problems  lacking  a  bilateral  symmetry, 
such  as  that  of  an  oblique  wing.  Furthermore,  the  small- 
disturbance  assumption  breaks  down  near  the  leading  edge; 
comparison  with  the  TSD  code  will  not  reveal  the  important 
limitations  of  the  present  (small-disturbance)  theory.  For  the 
above  two  reasons,  comparison  with  the  TSD  codes  has  not 
been  made.  Nevertheless,  we  believe  such  a  comparison  could 
be  quite  useful  in  future  studies  involving  shocks. 

From  a  theoretical  viewpoint,  a  limitation  of  current  3-D 
full-potential  codes  appears  to  be  the  empiricism  introduced 
by  the  modeling  of  the  inviscid  wake,  which  we  believe  could 
have  been  avoided.  Two  assumptions  on  the  vortex  sheet  were 
implicit  in  Refs.  25-28:  1)  the  shape  of  the  trailing  vortex  sheet 
is  specified  a  priori;  and  2)  a  condition  equivalent  to 


is  applied  at  the  TV  sheet.**  Assuming  a  thin  airfoil  section, 
assumption  1  is  seen  to  be  equivalent  to  that  in  the  TSD 
theory,  in  which  the  TV  sheet  is  transferred  to  the  “wing 
trace”  on  the .v axis.  Condition  2  is  incorrect  strictly  speaking, 
since  3-’0/3c2  is  generally  discontinuous  at  the  sheet.  In  fact, 
Eq.  (23)  would  lead  to 


far  downstream,  which  would  admit  only  span  loadings 
represented  by  a  second-degree  polynomial.  Apparently,  a 
central-difference  approximation  was  applied  to  Eq.  (23)  in 
Refs.  25-28,  stipulating  that  320/3c2  was  continuous  at  the 
sheet.  This  then  gives 


>T2 

/30X  - 

V dz  /  i\ 

V  dZ  /  TV 

that  is,  the  upwash  is  continuous  across  the  wing  trace,  which 
is  the  correct  requirement  according  to  the  TSD  theory. 
Therefore,  current  codes  based  on  Refs.  25-28  do  not 
correspond  to  the  exact  full-potential  solution. 

Compulations  of  i0,  o, .  and  e. 

The  reduced  mixed-type  problems  for  the  similarity 
solutions  0„,  0,,  and  0,  are  solved  numerically  by  a 
relaxation  method,  using  type-sensitive  difference  operators 
corresponding  to  Murman's  “fully  conservative  form."42 

1  There  is  a  similar  question  on  truncation  errors  near  the  symmetry 
plane  of  a  symmetric  swept  wing,  which  are  amplified  by  the 
singularity  of  the  coordinate  transformation. 

"This  controversial  condition  was  brought  to  our  attention  by 
Norman  1)  Malmuth. 
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The  procedure  of  line  relaxation  used  for  <i>0  may  be  con¬ 
sidered  standard,  except  for  the  uses  of  an  improved  far-field 
description  and  a  third-order  convergence  acceleration 
scheme. 45  Inasmuch  as  the  cases  studied  below  are  shock- 
free,  the  shock-fitting  algorithm44  and  shock -perturbation 
analysis  will  not  be  discussed.  A  more  complete  description  of 
the  solution  procedure  and  solution  examples  for  0O,  and 
<t>,  in  the  shock-free  cases  have  been  given  in  Ref.  16  (see  also 
Sec.  5.2  of  Ref.  34). 

Fsamples:  Comparison  with  Solutions  Based  on  the  Full-Potential 
PDF. 

For  the  purpose  of  comparing  the  lifting-line  analysis  (the 
similarity  solution  structure,  in  particular)  with  the 
corresponding  full-potential  solutions,  we  will  now  consider 
examples  of  high  subcritical  and  slightly  supercritical  com¬ 
ponent  (lows  over  oblique  as  well  as  symmetric  swept  wings. 
The  wing  planforms  considered  are  all  elliptic,  with  the  major 
axis  coinciding  with  the  midchord  (50%);  the  wing  sections 
are  generated  from  a  single  profile  NASA  3612-02,  40, 
rescaled  to  an  arbitrary  thickness  r;  the  latter  is  equated  to  a 
in  the  theory.  These  wings  fulfill  the  geometrical  description 
of  Eqs.  (8)  and  (19);  therefore,  the  inner  solution  can  be 
obtained  by  a  linear  combination  of  the  basic  similarity 
solutions  <j>,,  and  which  will  be  determined  after 
specifying  the  component  transonic  parameters  A'„,  the  ratio 
of  incidence  to  the  thickness  ratio,  and  the  locations  of  the 
leading  and  trailing  edges.  Note  that  this  airfoil  has  been  used 
in  various  wind  tunnels  and  preliminary  design  studies  of 
oblique  wings  at  =  0.60-1 .4. 11  Several  sets  of  basic 
similarity  solutions  have  been  obtained  for  A„  =  3.6  and  3.45 
for  this  airfoil,  and  have  been  described  in  some  detail  in 
Refs.  12-16.  These  solutions  are  used  in  the  subsequent 
comparison  with  the  full-potential  solutions.  The  case  with 
A„  =3.45  is  slightly  supercritical. 

Numerical  results  comparable  with  our  solutions  are 
generated  from  one  version  of  Jameson-Caughey  3-D,  full- 
potential  computer  codes,  “Fl.O  22”  (see  Refs.  6,  25,  and 
26),  which  was  modified  slightly  for  oblique-wing  analyses  at 


planform  K  a  14:1  ellipse:  the  wing  section  is  NASA  3612-02, 40,  eero 
incidence,  eero  twist,  and  eero  wing  Itend. 


NASA  Ames  Research  Center  and  Grumman  Aerospace 
Corporation.  The  different  algorithms  employed  in  FLO  22 
are  not  fully  conservative,  but  this  may  not  be  essential  for 
shock-free  solutions  presented  below.  We  point  out  that  the 
FLO  22  data  from  NASA  Ames  and  Grumman  are  not 
identical,  owing  to  the  differences  in  the  mesh  size,  spacing  of 
the  span  stations,  convergence  criteria,  and  t1'"  detail  of  the 
leading-edge  description.  The  availability  of  data  from  two 
sources  is  helpful  in  delineating  the  nature  of  discrepancy 
between  our  theory  and  the  more  exact  3-D  programs. 

A  number  of  FLO  22  runs  have  been  made  for  the  oblique 
wings  with  freestream  Mach  number,  swept  angle,  wing 
thickness,  etc.,  chosen  to  give  either  K„  =3.6  or  K„  =3.45, 
employing  the  same  basic  airfoil  section.  An  elliptic  planform 
is  used  in  each  case;  wing  twist  and  wing  bend  are  assumed  to 
be  zero.  Among  the  first  comparison  studies  made  is  a  case 
with  relatively  thin  wing  section  (6%  thickness  ratio)  and  a 
rather  high  aspect  ratio  (major-to-minor  axes  ratio  of  20).  The 
surface  pressure  coefficients  in  this  case  have  been  presented 
in  Ref.  15  (see  also  Fig.  9  of  Ref.  34).  The  freestream  Mach 
number  and  the  sweep  angle  used  in  the  FLO  22  calculation 
therein  are  =  0.8242  and  A  =  22.5  deg,  giving  a  component 
Mach  number  Mn  =0.7615.  Thus,  one  has  Kh  =3.60, 
0  =  1 .003,  and  t  =  0.1277  in  this  case;  the  component  flow  is 
slightly  below  being  critical.  There,  the  FLO  22  data  from 
NASA  Ames  and  Grumman  appear  to  be  rather  close  except 
next  to  the  leading  edge.  The  Cp  values  computed  from  the 
similarity  solutions  agree  reasonably  well  with  the  FLO  22 
data  except  near  the  leading  edge. 

The  degree  of  agreement  between  the  FLO  22  data  and  the 
theory  docs  not  appear  to  deteriorate  much  with  increasing 
wing  thickness  or  reducing  wing  aspect  ratio.  The  con¬ 
solidated  plot  in  Fig.  5  shows  the  surface  pressure  coefficient 
at  seven  span  stations  on  a  12%  thick,  14:1  elliptic  wing, 
pivoted  at  A  =  30  deg  for  =  0.7677.  The  component  Mach 
number  is  M„-  0.6648  in  this  case,  giving  A'„  =  3.45.  The 


Mich  number  of  0.7$$.  The  basic  plinform  Is  •  14:1  elHpse;  the  wing 
section  is  NASA  36124)2,  40,  rrro  incidence,  zero  (wisl,  ind  no  wing 
bend. 
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three  sets  of  surface  data,  ioi,  £;j,  and  d>_,(,  e  used  to 
construct  the  Cp  distributions.  The  Cp  value  for  the  com¬ 
ponent  flow  is  -0.689  and  both  FLO  22  and  lifting-line 
solutions  give  the  appearance  of  supercritical  shock-free 
regions  on  the  upper  surface  of  the  downstream  wing  panel. 

The  agreement  between  the  FLO  22  data  and  those  based  on 
the  asymptotic  analysis  in  Refs.  15  and  16  and  in  Fig.  5  should 
be  considered  as  being  better  than  expected,  inasmuch  as  the 
relative  error  in  the  asymptotic  theory  belongs  to  an  order 
determined  by  r-’-'  or  t whichever  is  larger.  The  magnitude 
of  t-'1  for  the  examples  shown  in  Fig.  5  is  0.243.  It  may  be 
recalled  in  this  connection  that  there  is  a  noticeable  difference 
between  the  two  sets  of  FLO  22  data  for  the  local  lift  coef¬ 
ficient  show  n  earlier  in  Fig.  4. 

It  is  not  altogether  clear  that  the  encouraging  agreement 
found  for  the  oblique  w  ings  may  still  hold  to  some  degree  with 
the  presence  of  an  apex  in  a  symmetric  swept  wing,  where  the 
theory  must  fail.  Figure  6  presents  a  consolidated  plot  for 
surface  pressure  coefficient  at  seven  span  stations  on  a 
symmetric  swept  wing  w  ith  the  same  basic  (elliptic)  planform, 
sweep  angle,  and  section  profile  as  those  in  the  preceding 
figure.  The  freestream  Mach  number  is,  however,  slightly 
lower  with  A/„  =  0.7549;  thus,  M„  =0.6538,  A'„  =  3.60, 
(  =0.1448,  and  0-1.062.  The  solution  remains  subcritical 
component-wise  for  most  stations  (the  dashed  line  in  Fig.  6 
indicates  Cr=  -  0.727,  w  here  the  local  component  Mach 
number  reaches  unity).  The  stations  v  =  0.92  and  0.97  are 
included  to  illustrate  the  solution  behavior  near  the  tip,  where 
the  trends  of  departure  of  the  asymptotic  (lifting  line)  from 
the  FLO  22  analyses,  as  well  as  between  the  two  sets  of  FLO 
22  data  themselves,  are  evident.  Away  from  the  tip, 
agreement  is  as  reasonable  as  in  the  two  preceding  com¬ 
parisons  for  the  oblique  wings.  However,  the  set  of  FLO  22 
data  in  filled  circles  gives  a  consistently  higher  peak  for  -C;, 
than  the  one  in  open  circles.  The  discrepancy  may  be 
associated  with  the  uncertainty  related  to  the  fluctuating 
spanwise  distributions  of  sectional  lift  shown  in  Fig.  4.  The 
present  theory  tends  to  give  a  consistently  lower  -  C(,  than  the 
FLO  data  on  the  lower  surface  around  the  quarter  chord.  This 
tendency  is  also  apparent  from  Fig.  5.  This  small,  but 
noticeable,  discrepancy  may  be  associated  with  the  accuracy 
and  degree  of  convergence  in  the  <t>„  solution. 

The  most  important  piece  of  information  from  Fig.  6  is  the 
comparison  made  for  the  station  closest  to  the  apex.  At  IOwo 
semispan  from  the  symmetry  plane  (.»’  =  0.097),  the  agreement 
of  the  asymptotic  results  with  the  FLO  22  data  (in  open 
circles)  remains  as  good  as  other  stations.  The  discrepancy  in 
the  vicinity  of  the  apex  is  seen  to  be  less  than  that  in  the  tip 
region:  the  latter  is  believed  to  arise  partly  from  the  treatment 
of  the  wing-tip  geometry  by  the  FLO  22  code  mentioned 
earlier.  In  passing,  we  remark  that  near  the  s  mmetry  plane 
where  the  leading-edge  sweep  angle  has  a  discontinuity,  even 
the  full-potential  code  may  not  be  strictly  reliable  owing  to  the 
scarcity  of  span  stations  (cf.  footnote,  third  paragraph  Sec. 
IV). 

V.  Concluding  Remarks 

The  foregoing  presentation  has  shown  that  Prandtl’s  lift¬ 
ing-line  idea,  originally  applied  to  an  unyawed,  straight  wing 
of  high  aspect  ratio,  can  be  extended  to  the  study  of  3-D 
mixed  flows  over  transonic  swept  w  ings.  Examples  of  oblique 
and  symmetric  swept  wings  involving  the  high  subcritical  and 
slightly  supercritical  component  flows  arc  shown:  com¬ 
parisons  with  full-potential  solutions  from  an  existing 
computer  code  are  made.  Except  near  the  wing  tip  and  root 
and  the  leading  edge,  where  breakdowns  in  either  (he  theory 
or  the  code  arc  expected,  the  agreement  with  the  full-potential 
results  is  encouraging. 

An  important  feature  of  the  analysis  is  the  upwash 
correction  proportional  to  )(df7dv),  resulting  from  the 
spanwise  0’’  )  component  of  the  locally  shedded  wake 

vorticity.  The  dependence  of  the  upwash  correction  on  the 


sweep  and  the  span  loading  have  been  shown.  A  significant 
gain  in  the  theoretical  development  through  this  approach  is 
the  availability  of  a  similitude  in  the  3-D  flow  structure  ap¬ 
plicable  to  wings  with  sections  generated  from  a  single  profile. 
With  the  similarity  solutions,  the  reduced  2-D  problems  are 
solved  only  once  for  all  span  stations.  From  these  and  the 
upwash  analysis,  the  dependence  of  the  flowfield  on  wing 
sweep,  aspect  ratio,  local  chord  variation,  as  well  as  the  twist 
and  the  bend,  can  be  explicitly  studied,  thus  furnishing  a 
rational  tool  for  transonic  aerodynamic  design. 

The  examples  studied  in  this  paper  do  not  include  examples 
with  embedded  shocks.  Solutions  making  use  of  the  3-D 
corrections  for  the  shock  jump  and  the  shock  geometry  have 
been  obtained,  and  their  comparison  with  full  potential 
results  will  be  discussed  in  a  sequel  paper.  Omitted  from  the 
foregoing  discussion  is  the  feasibility  of  solutions  via  a  special 
unsteady  analogy.  Unpublished  results  computed  from  this 
analogy  with  an  ADI  algorithm  by  T.  Evans  of  the  University 
of  East  Anglia  compare  rather  closely  with  similarity 
solutions,  and  promise  an  effective  procedure  for  treating 
more  complicated  shock  patterns. 

Whereas  local  treatment  of  the  leading-edge  breakdown  in 
the  small-disturbance  theory  is  possible  46  an  alternative 
solution  to  this  difficulty  is  to  develop  a  lifting-line  theory  on 
the  basis  of  the  potential  theory  without  the  small-disturbance 
assumption.  Work  along  this  line  has  been  carried  out  in  Ref. 
36,  where,  interestingly,  solutions  via  similarity  flow  structure 
may  again  be  found.  Treatment  of  nonuniformities  at  the  root 
and  along  the  tip  are  the  remaining  problems  for  the  lifting¬ 
line  theory  in  linear  and  nonlinear  regimes:  their  analyses  for 
the  transonic  regime  would  involve  solving  a  fuller  reduced  3- 
D  problem. 
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